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^ Abstract. A BCS model characterized by a phenomenological pair potential with 

on-site (Vq), nearest (V^), and next nearest (Vj) neighbour coupling constants, and an 
empirical quasiparticle dispersion taken from angle-resolved photoemission spectra 
is considered. The model can consistently explain the experimental data concerning 
^ the pair state of the hole doped cuprates. Three ingredients are required to make 

CN the interpretation possible: the existence of flat bands, a very small effective on- 

site repulsion, and a slightly dominating effective nnn attraction V2 of the order of 
§ 60-80meV with a ratio V2/V1 » 1.5. 

a 

X3 PACS Numbers: 74.20.-z, 74.20. Fg, 74.20. Mn, 74.72.-h 
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A number of recent experiments[l, 2] have indicated that the order parameter (OP) of the 
hole doped high-r^cuprates has a strong k dependence, and is very likely exhibiting line nodes. 
Moreover, a subset of these experiments [3], all of them using YBa2Cu307_^(Y123) samples, 
provided strong direct evidence that the OP changes sign under a 90° rotation, reminiscent of 
<i;f2_^2 -symmetry. However, another set of experiments on Y123[4] which also directly probe 
the k dependence of the OP show results that are hard to reconcile assuming a pure d^i^^ gap. 
Also, the temperature dependence of the London penetration depth measured on crystals of the 
electron doped Ndj g5Ceo j5Cu04_y(NCCO) is consistent with a weakly anisotropic BCS-like 
energy gap[5]. 

More recently, new information on the k dependence of the superconducting gap of high 
quality Bi2Sr2CaCu20g (Bi2212) crystals was reported from angle-resolved photoemission 
(ARPES) data by the Argonne group [6]. The gap on the Fermi surface (FS) was found to 
vanish at two points symmetrically displaced away from the (1,1) direction, but was finite on 
the (1,1) Une. This result is incompatible with a d-^i_yi-OV, which must have a single node 
along the (1,1) direction (a small admixture of an s-wave component, could slightly displace 
the node to one side, but cannot produce two nodes)[7]. Instead, the angular dependence of 
the gap resembled very closely that of an OP of the form Aj^ ~ cos cos ky [8], which we shall 
call the state. 

Considering these findings, one might draw the unpleasant conclusion that the pairing 
mechanism of the cuprates might be non-universal. A possible alternative was proposed by 
Chakravarty etal. [9] . In their theory, the main ingredient to the high is the postulate that the 
normal state of the cuprates is a Luttinger liquid confined to single Cu02-sheets. This implies 
non-coherent single-particle charge transport in the c-direction. In the superconducting state, 
however, coherent pair tunneling becomes possible, and the corresponding large gain in kinetic 
energy leads to the high T^.. Within this scenario, the symmetry of the pair state is determined 
by weak residual interactions within a single layer, and can thus easily vary between different 
materials. 

In this article, we make an attempt to find a consistent interpretation of the data within 
a simple BCS framework. We construct a phenomenological model that can qualitatively 
account for the apparently different pair states. The experimental quasiparticle dispersion 
relation, as determined by ARPES, serves as an input to the theory, and for the pair potential, 
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we allow on-site, nearest (nn), and next nearest neighbour (nnn) interactions. 

The main finding of the theory is the following: the degree to which the Cooper pairs can 
take advantage of the large single particle density of states (DOS) (arising from the observed 
flat bands in the quasiparticle dispersion of hole doped cuprates [10]), highly depends on the 
k dependence of the pair wavefunction. Assuming singlet pairing, and excluding the isotropic 
5-wave case by using a repulsive on-site interaction, we shall show that due to this effect, 
for the hole doped materials, there are precisely two different competing pair states, namely 
the d^2_y2 , and the s^^ state. Their relative stability is determined mainly by the ratio of the 
strengths for nn and V2 for nnn terms of the pair potential. A consistent interpretation of 
the experimental data can be achieved with a set of parameters for which the two states are 
almost degenerate. However, we find that the observed gap of electron doped materials can 
only be reproduced by an attractive on-site potential, which would indeed suggest that the 
pairing mechanism is different in this case. 

We consider the following model Hamiltonian on a square lattice 



cl-fj creating a quasiparticle of momentum k, spin a, where H^^^ is chosen to have the following 
form 



Here h stands for the number operator, d^id^J represents the nn (nnn) lattice vectors, and o is 
the spin index. The interaction constants Vi are free parameters, however, generally we assume 
Vq > 0, and Vp V2 < 0. This choice excludes isotropic s-wave pairing. 

The quasiparticle dispersion relation e(k) is obtained by a tight binding fit to normal state 
ARPES data on Bi2212, adjusted so that = (see Ref. [8] for details). Even though 
this fit was obtained from a particular material, we would like to stress that the qualitative 
features of the dispersion are common to all measured (by ARPES) high-r^materials near 
optimal doping. The conclusions of this paper do not depend on the detailed form of e(k), 
as long as the dominant feature is incorporated, the presence of flat bands in a large region 
of the Brillouin zone, leading to a van Hove singularity (vHS) in the DOS (see Fig. 1). The 
singularity seems to be sUghtly (= 30 meV for Bi2212, = 20 meV for Y123) below Ep in the 
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Table I: The basis functions 7],(k), their C^i, group representations, and the corresponding coupling constants entering 
the pair potential V(k - k') = J] to V', ?7Kk)?7/(k')- 



Ai 1 

1 A J ^(cosk^ + cosk^) 

2 5j i(cos^^-cos^p 

3 Aj cos^^cosfcy 

hole doped cuprates[10]. In the electron doped NCCO, flat bands are observed at = 300 meV 
below [11]. A fit of e(k) to the latter spectra looks similar to the one obtained for the hole 

doped materials (taking into account a shift in the chemical potential). 
After transforming to k space 

^int ~ 2A ^ ^^'^^ ^k+q,<j^^k'-q,CT'^k',C7'^^k,CT ' (3) 
q,k,k' 
a, a' 

where V(q) - Vq + (cos + cos q^) + SV2 cos q^ cos qy and A the volume, we apply standard 
BCS theory[12] to obtain the familiar gap equation 




Here E^- y £l + is the usual quasiparticle energy in the superconducting state, /3 the inverse 
temperature. The potential V(k - k') is separable, and restricting ourselves to singlet pairing, 
can be written as V{k - k') = ^^^q V,T]j(k)T],(k'). Each basis function T]j(k) can be associated 
with one of the irreducible representations (in this case, all of them one-dimensional) of the 
C^^ point group which characterizes the square lattice. They are listed in table II. The OP is 
expanded as = ^to ^i^^)- 

As usual, is determined by linearizing Eq. (4), which we can then solve separately 
within each irreducible representation. The representation with the largest T^. is the one which 
will be chosen by the system, and as long as there is no additional phase transition, the OP at 
r < is given by the continous evolution of this solution as T is lowered. Here, we shall not 
consider the consequences of a second transition. 
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The linearized gap equation becomes one-dimensional for the representations (d^2_y2) 
and B2 (4ry)» however for the case, we have mixing of the three basis functions T7,(k), i = 
0, 1, 3, and thus three coupled equations. In each irreducible representation fH, the linearized 
gap equation can be written as 

At this point it is instructive to define a new quantity for each basis function rij(k) as = 
X ^^(k)5(e-ek), which we call its pairing density of states (PDOS). The role of the PDOS in 
the gap equation for an anisotiopic pair state is analogous to the role of the usual single particle 
DOS, A^(e), for an isoti-opic OP (for T7o(k) = 1, P^ie) = Nie)). The PDOS of a particular pair 
state is of crucial importance, since it is the value of Pf{e) around e^, which (apart from the 
pair potential strength) determines its stabiUty. 

In Fig. 1, we plot P,(e) for the basis functions under consideration. This plot reveals an 
interesting and important point. The vHS in the DOS is reflected in the PDOS of only two 
of the four anisotropic basis functions, namely the d^2_y2, and the s^ state. The s* state has 
almost vanishing weight, and the d^ state is also not very strong at e^. From this it is already 
clear that there are two competing pair states in the relevant doping range of the hole-doped 
cuprates. Their relative stabiUty is determined mainly by the ratio V1/V2. 

Now let us discuss in more detail the solutions of Eq. (5). Our strategy is as follows: first, 
we determine a set of realistic parameters for the coupling constants which can reproduce the 
measured k dependence of the gap in Bi2212, i. e., an almost pure gap. Second, we see, 
whether the same or slightly modified set of parameters gives results that are consistent with 
data from the other cuprates. We solve Eq. (5) numerically on a dxd k-space grid, where 
the number d at which convergence is obtained varies between 500 and 8000 depending on 
the size of the coupling constant. If the coupUng constant is small, so is T^, and one has to 
use more k-points to accurately sample the important region around the FS, where the factor 
e^^ tanh(j8eij/2) is large. Note that the usual cutoff that prevents the logarithmic divergence in 
standard BCS theory is not necessary here due to the finite bandwidth. 

The first criterion necessary to reproduce the Bi2212 gap is that an solution be stable 
versus the d^ state. Since the PDOS of the s* state is so small around e^, its component in the 
OP is negUgible, unless we consider values IVjl > IVjL for which the d^2_y2 state is clearly the 
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Figure 1 : The PDOS P(e) for the various OPs considered. In case of the isotropic gap P(e) = N{e), the 
DOS. The inset shows the experimental Fermi surfaces for Bi2212, and NCCO as obtained from a tight 
binding fit to ARPES data. The wide region between the dotted and the continous fine encloses the flat 
bands around the vHS (Evhs = -34meV). 

most favourable. Hence, in this comparison, we can safely set = 0. To gain some insight, 
we investigate the phase boundary between the two states in the Vq, V2 plane. Figs. 2(a-c) 
show the phase boundaries for hole doping concentrations 5 = 0.17 - 0.26 (note that the vHS 
corresponds to 5 = 0.34). At Vq - 0, we find that the Aj state is favoured for IV2I ^ 200meV. 
The stability of the d^y state at large coupling is explained as follows. The effective region 
being sampled in the equation grows as increases [this is easily understood from the jS^ 
dependence of the factor e^^ tanh(P^eJ2)]. Hence, the contribution to the pairing by the 
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Figure 2 : (a-c) The phase boundary between the dxy and the most stable Ai solution for three different 
doping levels (Vi = 0). We did not determine the boundary, in the white region of small IV2I, since 
Tc < IK. (d) The angular dependence of the T = superconducting gap in the Ai representation for 
various strengths of the on-site repulsion Vq, calculated for the experimental FS of Bi2212. The coupling 
constants (in meV) were chosen as Vi = 0, (Vq, V2) = [(0, -60), (100, -72), (200, -80), (400, -90)], which 
results in a value Tc = lOOK for all four sets. 

region around e = 0.15, where its PDOS is large, grows substantially with increasing IVjI- 
For Vq > 0, the region of Aj stability narrows to a IVjl stripe centered at about lOOmeV. As 
a function of increased hole doping, the phase boundary moves substantially towards larger 
Vq. This is a result of the increasing ratio of the to the PDOS [/'3(e^)//'4(e^)] as we 

approach the vHS (see Fig. 1). The region of realistic values for < 200K is obtained for 
IV2I ^ lOOmeV regardless of Vq. This means that the Aj solution is indeed stable versus the 
state in the physically interesting parameter regime. 

The second criterion that needs to be satisfied is the k dependence of the gap, in partic- 
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ular, the position of its nodes, and its relative magnitude on either side of the nodes. Within 
the Aj representation, this is determined by the ratio Vq/V2. On general grounds, one can al- 
ready anticipate the result: Since the penalty of a large Coulumb on-site repulsion can only 
be avoided if the two electrons making up a Cooper pair never visit the same site at the same 
time, one has to require that Y^ii'^^'^u) - 0- In BCS theory, this condition is equivalent to 
^^^(r-r' = 0,t = 0) ~ Xlk^k^^k - (^^ the anomalous propagator), and requires a sign 
change of the gap, as well as roughly equal positive and negative portions of it in the region 
near the FS which contributes strongest to the sum. Hence for large Vq, we expect the nodes to 
converge to an angle near = 7t/S (22.5°). A quantitative analysis confirms this argument. In 
Fig. 3(d), we plot the (p dependence (measured from the [k, k] point) of the superconducting 
gap at doping 5 = 0.17 (the experimental Bi2212 value) on the FS for various values of the 
on-site repulsion Vq. For each Vq, V2 was chosen, so that ~ lOOK. For a value as low as 
Vq = lOOmeV, we already find a shift of the node by 5°, and the size of the gap at 45° has more 
than doubled. Consequently, we have to conclude that the effective on-site repulsion between 
the quasi-particles has to be very weak to explain the observed gap in Bi2212 within an Aj 
scenario. This is certainly a surprising result which requires a microscopic justification, since 
the bare copper on-site repulsion is generally believed to be large. 

In the remainder of this paper, we concentrate on the competition between the and the 
d^2_y2 State. As a consequence of the previous considerations, we set VQ = Om the calculations 
that follow. Fig. 3(a) shows the phase boundary between the two OPs in the relevant doping 
region. The boundary is essentially linear, and defined by a ratio of V2/V1 ~ 1.4. As a function 
of increasing hole doping there is a slight shift of the line towards smaller values of 1^21- This 
means that in principle a transition from a d^2_y2 to a state is possible when the doping 
is increased provided that the two solutions are almost degenerate. Note, that this is again 
a straightforward consequence of the sUghtly different e dependence of the accompanying 
PDOS. To illustrate this point further, in Fig. 3(b) we plot the doping dependence of for the 
two states using two different sets of parameters. We note, that if the state is stable, it is 
favoured only in the vicinity of the vHS. A transition to the d^2_y2 state has to occur on either 
side of the optimal doping regime, provided \Vi\ is not extremely small. Figs. 3(c,d) show 
the coupUng dependence of for the two competing states at two different doping levels 5 = 
0. 17, 0.26. At large coupUng the dependence becomes approximately Unear in Vp V2 resp.[13], 
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Figure 3 : (a) The phase boundary between the dx^-y^ and the most stable A i solution for three different 
doping levels (Vq = 0). (b) The doping dependence of for the dx^-y^ and the A\ state calculated using 
two different sets of parameters, p\,p2 given by (Vi, Vj) = [(-50, -70), (-40, -58)] (in meV). (c,d) The 
coupling dependence of Tc for the two states at two different doping levels 5 = (0.17, 0.26). 

a consequence of the proximity to the vHS. For > 20 - 30meV the integrated PDOS over 

the region close to £p, where e|^^ tanh(/3^e,j/2) ~ /3^/2, essentially saturates as a function of 

increasing j8^, since it is exhausted by the vHS (see Fig. 1). Hence the dependence of the 

sum in Eq. (5) becomes approximately linear. 

Until now, we were mainly concerned about the determination of a set of parameters which 

can reproduce the observed gap in Bi2212. In the final part of the paper, we want to see 

whether a similar set of parameters can consistently explain the experimental data relevant to 

the OP symmetry in Y 123. As already noted by other authors [14], the lower orthorhombic 

symmetry of Y 123 has the effect that the five basis functions we considered now belong to 

only two different irreducible representations of the relevant point group C2„. In particular, 
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the d^2_y2 state now falls into the identity representation Aj, and a mixing with the s-wave like 
basis functions at is now possible. Note, that even though the true symmetry of Bi2212 is 
also orthorhombic due to the superlattice modulation, since the latter occurs along the (1,1) 
direction, the <i^2_^,2 state does not belong to the representation in this case, and hence no 
mixing is possible[8]. To simulate the orthorhombic distortion, we artificially introduce a 
difference in the versus ky quasiparticle dispersion characterized by a parameter a defined 
by setting ,^ = (1 - = (1 + a)?,. Here is the real space hopping matrix element for 

a lattice vector 5, in the tight-binding formulation. The quaUtative effect of this symmetry 
breaking does not depend much on the value of a, so we fix it at a = 0.05. In our model, the 
orthorhombic distortion splits the vHS into two peaks, roughly separated by 50meV. The peak 
at higher energy corresponds to doping d = 0.24. Note, that this splitting is of no relevance to 
the arguments presented below. 

Since we already learned that the two competing states in the hole-doped region are the 
d^2_y2 and the state, in the following analysis, we restrict ourselves to the mixing of those 
two states, i. e., we write the gap function as A^^(T} - Ai(r)r]2(k) + AjiTjrj^Qi.), and now 
we solve the full non-linear gap equation to determine the coefficents A]^(T), A2(T). In the 
parameter regime previously determined (where d^2_y2 and s^y are almost degenerate) , we find 
considerable mixing between the two states. In Figs. 4(a,b), we show the angular dependence 
of the r = OP at two different doping levels 5 = 0.24, 0. 19 for various values of the coupUng 
constants Vp V2 going from the case of a dominating to a dominating d^_y2 component. 
The crucial result here is that at T = 0, the two solutions repel each other, i. e., there is no value 
of the coupling constants where Aj = The maximum admixture of the minor component is 
roughly 25%. This is clearly demonstrated in Fig. 4(d), where the temperature dependence of 
the two components is plotted for two sets of parameters. For T close to T^, the two components 
develop roughly equally with decreasing T until at some T, one of the components starts to be 
boosted, and the other one becomes suppressed. This cross-over temperature can be related to 
the crossing of the free energy of the two solutions in tetragonal symmetry, which can occur 
when they have a similar T^. If the solutions are in the same representation, however, a crossing 
cannot occur; the solutions come very close to each other, and at the cross-over temperature, 
the stable solution changes its character rapidly, but continously. 

This behaviour can lead to an interesting effect, which is demonstrated in Fig. 4(c). It is 
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Figure 4 : Effect of the orthorhombic distortion: (a,b) The angular dependence of the T = mixed 

d;^-^ + Sxy OP in orthorhombic symmetry at two doping levels for different sets of parameters, 5 = 0.24, 
/?i • • -PS defined by Vi = -40, = [-65, -60, -57.6, -57.2, -40], and 5 = 0. 19, pi---p5 defined by 
Vi = -40, V2 = [-65, -60, -58.8, -58.5, -40]. (c) The angular dependence of the mixed d.i.yi + s„ OP at 
different temperatures (t = TITc). (d) The temperature dependence of the d^^-y^ component Ai, and the Sxy 
component A2 of the mixed d^^-yi + s^y OP for two different sets of parameters, p\, p2 given by Vi = -40, 
y2 = [-57.5,-57]. 

actually possible that as a function of T, the solution changes so dramatically, that close to T^, 
the sign of the gap is the same at 0° and 90°, whereas as T is lowered, the d^2_y2 character 
starts to dominate, and the gap changes sign under a 90° rotation. Experimental evidence for 
such behaviour has already been reported on samples of Y 123 [15]. 

In any case, for the set of parameters determined to fit the gap of Bi2212, we find consider- 
able mixing of the d^2_y2 and the state in Y 123. As shown, it is not difficult to obtain a gap 
function with predominantly d^2_y2 character so that the OP changes sign under a 90° rotation, 
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while its average over the FS is still considerably large due to the component. This could 
resolve the controversy which arose from results of different phase sensitive measurements on 
Y123 mentioned at the beginning. 

In conclusion, we have shown that a phenomenological model based on standard BCS 
theory can consistently explain the experimental data concerning the pair state of the hole 
doped cuprates. Three ingredients are required to make this interpretation possible: (i) The 
existence of flat bands (ii) a very small effective Vq (iii) a slightly dominating effective nnn 
attraction of the order of 60-80meV with a ratio Vj/V^ = 1.3-1.5. Clearly a derivation of 
the effective model presented here from a Hamiltonian with reaUstic, physically motivated 
parameters and interactions is an important next step in the understanding of the true pairing 
mechanism. Finally, we just briefly mention that the electron doped NCCO does not fit into 
this scenario. Given the experimental quasiparticle dispersion we found that only isotropic 
s-wave pairing, i. e., a negative V,, = -200meV, can reproduce the experimentally observed 
nodeless superconductivity. For the anisotropic states considered, the d^i^yi state is by far the 
most stable in this case. 

We acknowledge financial support by the National Science Foundation (NSF-DMR-91- 
20000) through the Science and Technology Center for Superconductivity (R.F.), and of the 
U.S. Department of Energy, Office of Basic Energy Sciences, under Contract No. W-31-109- 
ENG-38 (M.R.N.). 
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